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Abstract 

We address the problem of coupling non-Hermitian systems, treated 
as fundamental rather than effective theories, to the electromagnetic 
field. In such theories the observables are not the x and p appearing 
in the Hamiltonian, but quantities X and P constructed by means of 
the metric operator. Following the analogous procedure of gauging a 
global symmetry in Hermitian quantum mechanics we find that the 
corresponding gauge transformation in X implies minimal substitu- 
tion in the form P — > P — eA(X). We discuss how the relevant matrix 
elements governing electromagnetic transitions may be calculated in 
the special case of the Swanson Hamiltonian, where the equivalent 
Hermitian Hamiltonian h is local, and in the more generic example of 
the imaginary cubic interaction, where H is local but h is not. 

PACS numbers: 03.65. Ca, ll.30.Er, 02.30.Mv 



1 Introduction 

Recent interest in Hamiltonians that are non-Hermitian but nonetheless have 
a real spectrum dates from the pioneering paper of Bender and BoettcherpQ, 
which gave strong numerical and analytical evidence that the spectrum of 
the class of Hamiltonians 

H = p 2 +mV - (ix) N (1) 

was compieteiy real and positive for N > 2, and attributed this reality to 
the (unbroken) PT symmetry of the Hamiltonian. Subsequently a large 
number of PT-symmetric models were explored (see, e.g. ¡2]), and it was 
found that the phenomenon was rather general. The drawback that the 
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natural metric on the Hilbert space, with overlap J ipi(—x)ipj(x)dx, was not 
positive definite, was overeóme by the realization[3] that one could construct 
an alternative, positive-definite metric involving the so-called C operator. 
The formalism was further developed by Mostafazadeh|4], building on earlier 
work by Scholtz et al. [5]. In particular he showedjü] that such a Hamiltonian 
H was related by a similarity transformation to an equivalent Hermitian 
Hamiltonian h. The key relation is the quasi-Hermiticity of H: 

H* = t]Ht}- x , (2) 

where r¡ is Hermitian and positive definite. r¡ is related to the C operator by 
rj = CP, and it is frequently extremely useful[7] to write it in the exponential 
form r¡ = e~® . Occasionally r¡ can be constructed exactly (see, for example ¡Hl 
El Uní CEU H2¡)> but more typically it can only be constructed in perturbation 
theory, for example for the ix 3 model[T3j. 

From Eq. (¡2]) we can immediately deduce that 

h = pHp~\ (3) 

is Hermitian, where p = e~^®. Other operators A will also be observables, 
having real eigenvalues, if they are also quasi-Hermitian, i.e. 

A ] = r]Ar]-\ (4) 

and they again are related by the similarity transformation to Hermitian 
counterparts a: 

A = p~ x ap . (5) 

The similarity transformation also transforms the states of the Hermitian 
system, \<p), to those of the quasi-Hermitian system, \ip): 

W)=p- X \<P). (6) 
This implies that the matrix element of an operator is 

(o)«,- = Mto0K->- (7) 

In particular, the matrix elements of an observable can be written as 

= (fM^j) ■ (8) 
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A very important observation is that 



H(x,p) 



H{pXp-\pPp- 1 ) 

P H(X,P)p- 1 

h(X,P). 



(9) 



Thus an alternative way of finding h is to calcúlate the observables X and 
P and then rewrite H(x,p) in terms of them. 

The above concerns quasi-Hermitian systems considered in isolation. How- 
ever, important conceptual issues arise when one attempts to consider such 
systems in interaction with an otherwise Hermitian environment. For ex- 
ample, Ref. [H] examined a non-Hermitian analogue of the Stern-Gerlach 
experiment in which the role of the intermedíate inhomogeneous magnetic 
field flipping the spin is taken over by an apparatus described by a non- 
Hermitian Hamiltonian. This type of set-up has been further discussed and 
elaborated in a series of papers by various authors[I5l HS1 El CHl EB l2"ü] . 

Again, scattering gives rise to problems, since unitarity, as conventionally 
defined, is generically not satisfied for a PT-symmetric Hamiltonian. Uni- 
tarity can be restored, by use of the r\ metric, but then the concept of "in" 
and "out" states has to be drastically Ref. ¡21], [22], or in some cases¡23] less 
drastically, revised. 

The present paper is concerned with another such issue, namely how one 
couples a charged particle described by a quasi-Hermitian Hamiltonian to the 
electromagnetic field, following as closely as possible the well-known gauging 
procedure for a Hermitian Hamiltonian. This problem has been previously 
dealt with by Fariah and Fring[23] in a treatment which in many ways is 
more sophisticated than the present paper, dealing with pulses rather than 
plañe waves and going beyond first-order perturbation theory. However, 
the subtleties arising from the difference between x and X (see Eq. ( 1221 
below) were not encountered there because the calculations were done entirely 
within the framework of the dipole approximation, where the electromagnetic 
potential A is just a function of time. 

2 Brief review of the standard procedure 

In standard quantum mechanics the probability density is just |?/>(íe)| 2 , which 
is unchanged under a change of phase of the wavefunction: í¡j — ► e iea ip pro- 
vided that a is a real constant. If we try to extend this to a = a(x), a 
real function of a?, an extra term appears in the Schródinger equation, be- 
cause now pe íea ip = e iea (p + eVa)ip. We cancel this additional Va term by 
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minimal substitution: 



p —> p - eA . (10) 

Then under the combined transformations 

i¡) — > tp' = e iea íjj , , 

A^A' = A-Va, 

we obtain (p — eA)ip — > e iea (p — eA)ip, as required. Moreover the electric 
and magnetic fields are unchanged by the gauge transformation (1221) . 
So for a normal Hamiltonian of the form 

H=^ + V(-), (12) 

the coupling to the vector potential is —e(A.p + p.A)/(2m). In first-order 
perturbation theory a standard procedure then gives the transition rate be- 
tween the states |¿) and \j) induced by a plañe wave 

A{x, t) = J duA{uj)e l{k - x - ut) + ce. (13) 

as 

e 2 

^•oc— \(i\ PA \j)\ 2 , (14) 

in the dipole approximation e lk ' x ~ 1 over the range of the interaction. Here 
the constant of proportionality is (2n/h 2 )A(uJij) 2 , where uj^ = — Ej)/h, 
and pa is the projection of p in the direction of A. 

The matrix element (i\pA¡j) can be recast in terms of {i\xA¡j), where xa 
is similarly defined, by 

(Ei - Ej){i\x\j) = (i\[H,x}\j) = --{i\p\j), (15) 

TXXi 

so that 

(í\Pa \j) = imujijiilxA \j) ■ (16) 



3 Quasi-Hermitian quantum mechanics 

The totafl probability is now (^l^l^), where t] is the metric operator. This is 
no longer invariant under \ip) — > e tea ^ except in the special case where 
r\ = T)(x), so that [77, x] = 0. 

TNote that the probability density g(x) = (tp\p\x) (x\p\i¡)) is also invariant under the 
transformation of Eq. (JTTJ) . 
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It is, however, invariant under 

W-^e^W, (17) 
where X is the observable X = p~ l xp. For then 

{ip\r]\*P) -> (^|e- íea(x)t r/e íea(x) |^) 

= (^%|^>, (18) 

since = Note that, in terms of the eigenstates \ip) of h, Eq. ( ÍT71) 
corresponds to 

|y?> -> pe ieQ(x) p"V) = e iea{x) \ip). (19) 

Since we are using X in the exponent in Eq. (j!7p . we will also need 
to write in terms of X and the corresponding conjúgate observable P, 
according to Eq. (¡HJ), i.e. 

H(x,p)=h(X,P). (20) 
The minimal substitution we require, in h(X,P), is then 

P^P-eA(X) (21) 

with the combined transformations 

^) _> |^') = e *«(-X)|^) 

A(X) -> A'(X) = A(X) - V x a(X) . 



'22) 



It is important to note that because X and x do not commute, the argument 
of A in Eq. (j2*T|) must be X rather than a; in order to ensure that 

e -iea(X)(p _ eA lyea(X) = p _ eA 

Given the gauge transformation of Eq. f)22p . we are obliged to define 
B(X) = Vx x A(X), and the Fourier transform of Eq. (ÍTBl will also have 
to be rewritten in terms of X. How are we to interpret this, when X is 
a complicated non-local operator? The answer is that the external, classi- 
cal electromagnetic potential is in reality A(£), where £ is a real vector of 
position. Then B(£) = x A(£), and Eq. ÍÍT5]) becomes 



A(£,í) = J duÁ{u)é ik ^- ut) + ce. 



(23) 



Then, in the interaction with the non-Hermitian system, £ is replaced by 
the operator X, of which it is the eigenvalue. This is in parallel with the 
normal practice whereby in Eq. ( ÍT3l) it is understood that a? is a numerical 
vector, but in its interaction with a Hermitian system x is interpreted as the 
operator x. 

If h is of standard form, p 2 /(2/i) + U(x), the scattering rate is 

e 2 

w i:j oc — \{i¡Ji\r¡P A \^3)\ 
/i 

= - 2 \(^\Pa\^}\ 2 , (24) 
/i 

and the second form of the matrix element can then be rewritten, as in the 
Hermitian matrix element of xa, namely 

(ipi\pA\ifj) =ÍflWij{tPí\xA\tPj)- (25) 



3.1 The Swanson model 

A much-studied example where h, but not H, is of standard form is the 
Swanson HamiltonianjB], whose three-dimensional versión reads 

p 2 1 1 

H = ^— + -iue{x r ,p r } + -m 2 uj 2 x 2 , (26) 

with m 2 = (1 — £ 2 )m\. There is actually a one-parameter family[2Tj of Qs, 
from which we consider just the two cases (i) Q = Q(x) and (ii) Q = Q(p). In 
either case the equivalent Hermitian Hamiltonian is just a harmonic oscillator 
of the form 

h(x,p) = ^- + ±fjw 2 x 2 (27) 

(i) Q = Q(x) = e m\uj x 2 . 

This amounts to completing the square as 

(p + iem-iux) 2 1 , . 

H = ^— '— + -miufx 2 28 

2mi 2 1 y ! 

so that X — x, while P = p + ie m\ux. Thus in this case 

p 2 1 

h(x, p) = J— + - mi cü 2 x 2 , (29) 
¿mi 2 



so that /i = mi. The coupling to the vector potential is thus 

A.P + P.A) = — [(A.p + p.A) + iem l uü{A.x + x.A)} (30) 



2m 1 2mi 
The required matrix element 

(^PaI^j) = ((pi\pA\<Pj) t (31) 

is then found from expressing each component of p on the right-hand side in 
terms of creation and annihilation operators: p = i\/(m\0J /2)(a' — a). 

(ii) Q = Q(p) = -e x 2 1 \m 2 uj) 

This amounts to completing the square instead as 



H = 1 — m 2 í¿; x + 



2m 2 2 \ m,2LJ 
1 



(32) 



2m 2 2 

so that P = p, while X = a; + zep/(m 2 u;). Thus in this case 

fc(sc,p) = 7^— + \m 2 uj 2 x 2 , (33) 
zm 2 2 

with fi = m 2 . The coupling to the vector potential is thus 

[A.P + P.A) = {A.p + p.A) (34) 



2m 2 2m 2 

The matrix elements are still of the form of Eq. (131 j) . but now the components 
of p on the right-hand side are expressed as p = i-\/(m 2 a;/2)(a^ — a). 

The important thing to notice is that one will get different transition 
rates in the two cases. That is, the system is determined not only by the 
Hamiltonian H, but also by the particular metric operator rj used to restore 
unitarity. 

3.2 Imaginary cubic interaction 

The more common situation is that H is of standard form, while h is a com- 
plicated non-local object. For example, in the case of the (one-dimensional) 
prototype Hamiltonian 

H = ^(p 2 + x 2 )+igx 3 , (35) 
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we have^ 

'4 



Q = -g[-p 3 + 2xpx)+0(g 3 ), (36) 



which gives rise¡26¡ [27j to the observables 

X = x + ig(x 2 + 2p 2 ) + g 2 (—x 3 + 2pxp) 

P = p — ig[xp + px) + g 2 (2p 3 — xpx) 



+ 0(g 3 ). (37) 



Referring to Eq. (~"~J), we can write H(x,p) as h(X, P), where h(x,p) has been 
calculated up to second order in g as ¡26], 



h(x,p) = ±(p* + x 2 ) + 3g 2 (±x* + S 2i2 (x,p)-~j+0(g 4 ), (38) 

where 82,2(2, p) = (x 2 p 2 + xp 2 x + p 2 x 2 )/3. 

From Eq. (¡3*51) . we see that the minimal substitution P — > P — eA(X) in 
h(X,P) will give rise to additional couplings, of order g 2 , arising from the 
mixed term S 2 ,2(X, P). 

To 0(g) the matrix elements will be just (^'qPj^ijjj). In order to calcú- 
late this we will need the corrected eigenfunctions, which have a first-order 
contribution, namely 

^(x) =^(x)+gJ2(^ 3 \^)^(^ + 0(g 2 ) (39) 
In this case it is much easier[2H] to work with H directly rather than with h. 



4 Summary 

For a standard Hermitian system the coupling to the electromagnetic poten- 
tial, via the minimal substitution p — > p—eA(x), is induced by implementing 
the position-dependent phase change ip — > e tea ^ip and demanding that the 
transformed Schrodinger equation be unchanged. For a quasi-Hermitian sys- 
tem we find instead that the phase must be taken as a(X), where X is the 
observable associated with x. The coupling to the electromagnetic vector 
potential thus induced is via the minimal substitution P — ► P — eA(X) in 
H(x,p) written in terms of X and P, where P is the observable associated 
with p. 

The matrix elements governing electromagnetic transitions from one state 
of the system to another depend on both H and the metric r¡. In the special 



case of the Swanson Hamiltonian, when the equivalent Hermitian Hamilto- 
nian h is local, this dependence is encoded in the mass of the particle, which 
cannot simply be read off from H. Generically h is not local, and the cou- 
pling is considerably more complicated. 

Acknowledgement 

I am grateful to the referee for pointing out a serious mistake in the first 
versión of this paper. 

References 

[I] Bender C M and Boettcher S 1998 Phys. Rev. Lett. 80 5243 
[2] Levai G and Znojil M 2000 J. Phys. A 33 7165 

[3] Bender C M, Brody D C and Jones H F 2002 Phys. Rev. Lett. 89 270401; 
ibid. 92 119902 (E) 

[4] Mostafazadeh A 2002 J. Math. Phys. 43 205 

[5] Scholtz F G, Geyer H B and Hahne F 1992 Ann. Phys. 213 74 

[6] Mostafazadeh A 2003 J. Phys. A 36 7081 

[7] Bender C N, Meisinger P N and Wang Q 2003 J. Phys. A 36 1973 
[8] Swanson M S 2004 J. Math. Phys. 45 585 

[9] Geyer H B, Scholtz F G and Snyman I 2004 Czech. J. Phys. 54 1069 
[10] Jones H F 2005 J. Phys. A 38 1741 

[II] Jones H F and Mateo J 2006 Phys. Rev. D 73 085002 
[12] Assis P and Fring A 2009 J. Phys. A 42 015203 

[13] Bender C M, Brody D C and Jones H F 2004 Phys. Rev. D 70 025001; 
ibid. 71 049901 (E) 

[14] Bender C M, Brody D C, Jones H F and Meister B K 2007 
Phys. Rev. Lett. 98 040403 

[15] Assis P and Fring A 2007 Preprint [arXiv:quant-ph/0703254| 

[16] Guenther U, Rotter I and Samsonov B F 2007 J. Phys. A 40 8815 

9 



[17] Mostafazadeh A 2007 Phys. Rev. Lett. 99 130502 

[18] Bender C M, Brody D C and Jones H F 2008 Preprint la7Xiv:0804.3487l 
[quant-ph] 

[19] Nesterov A 2008 Preprint larXiv:0806.4646l [quant-ph] 

[20] Guenther U and Samsonov B F 2008 Preprmt laj5Eiv:0807.3643V 2 [quant- 
ph] 

[21] Jones H F 2007 Phys. Rev. D 76 125003 

[22] Jones H F 2008 Phys. Rev. D 78 065032 

[23] Znojil M 2008 Preprint larXiv:0805.2800l [hep-th] 

[24] Faria C and Fring A 2007 Láser Physics 17 424. 

[25] Musumbu D P, Geyer H B and Heiss W D 2007 J. Phys. A 40 F75 

[26] Jones H F 2005 J. Phys. A 38 1741 

[27] Mostafazadeh A 2005 J. Phys. A 38 6557; ibid. 8185 (E) 

[28] Bender C M, Chen J-H and Milton K A 2006 J. Phys. A 39 1657 



10 



